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ABSTRACT 

We present a highly complete and reliable mid-infrared (MIR) colour selection of 
luminous AGN candidates using the 3.4, 4.6, and 12 (im bands of the WISE survey. 
The MIR colour wedge was defined using the wide-angle Bright Ultra-Hard XMM- 
Newton Survey (BUXS), one of the largest complete flux-limited samples of bright 
(/4.5-10 kev > exlCT^ergs^cnr 2 ) "ultra-hard" (4.5-10 keV) X-ray selected AGN to 
date. BUXS includes 258 objects detected over a total sky area of 44.43 deg 2 of which 
251 are spectroscopically identified and classified, with 145 being type-1 AGN and 
106 type-2 AGN. Our technique is designed to select objects with red MIR power- 
law spectral energy distributions (SED) in the three shortest bands of WISE and 
properly accounts for the errors in the photometry and deviations of the MIR SEDs 
from a pure power-law. The completeness of the MIR selection is a strong function 
of luminosity. At L2-iokev > 10 44 ergs -1 , where the AGN is expected to dominate 
the MIR emission, 97.1±H% and 76.5^111% of the BUXS type-1 and type-2 AGN 
meet the selection. Our technique shows one of the highest reliability and efficiency 
of detection of the X-ray selected luminous AGN population with WISE amongst 
those in the literature. In the area covered by the BUXS survey our selection identifies 
2755 AGN candidates detected with SNR^5 in the three shorter wavelength bands of 
WISE with 38.5% having a detection at 2-10 keV X-ray energies. We also analyzed the 
possibility of including the 22/j.m WISE band to select AGN candidates, but neither 
the completeness nor the reliability of the selection improves. This is likely due to 
both the significantly shallower depth at 22/mi compared with the first three bands 
of WISE and star-formation contributing to the 22/zm emission at the WISE 22/im 
sensitivity. 
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1 INTRODUCTION 

There is strong observational evidence that active galactic 
nuclei (AGN) play an important role in the f ormation and 
growth of galaxies (e.g. iMagorrian et al.lll998l ). Most super- 
massive black hole growth takes place during an obscured 
quasar phase, as suggested by the integrated energy density 
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of the cosmic X-ray background jFabian fc Iwasawal [19991 ). 
To understand the evolution of galaxies and to trace the en- 
ergy output due to accretion and its cosmological evolution, 
it is critical to map the history of obscured accretion. 

X-ray surveys with XMM-Newton and Chandra at 
energies <10 keV are sensi tive to all but the most 
heavily obscured AGN (e.g. iDella Ceca et all 120081 ). In 
Compton-thick AGN (rest-frame column densities exceed- 
ing N H ~ 1.5 x 10 24 cm~ 2 ) the observed flux below 10 keV 
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can be as low as a few % of the intrinsic nuclear flux. In the 
Compton-thick regime the high energy photons that survive 
the photoelectric absorption get scattered in the absorber 
losing part of their energy (Compton down-scattering). 
This is an important effect that can significantly suppress 
the transmitted continuum (|Mattll2002l ; (Murphy fc Yaaoobl 
120091 ; lYaaoob et al] |2010| ). The ongoing Swift/BAT and 
INTEGRAL/IBIS all-sky surveys at energies 15-200 keV 
are providing the least bia sed samples of absorbed AGN 



th e local Universe (e.g. Bird et all 20071 ; iTueller et al.l 



120081 ; IWinter et aT]|2009l ; iBurlon et a.l.ll201lh . However 



these surveys are bias ed against the mos t heavily absorbed 
Compton-thick AGN |Burlon et alj|201lh . 

Surveys at mid-infrared (hereafter MIR) wavelengths 
(>5/im) are much less affected by extinction since the ob- 
scuring dust re-emits the nuclear optical-to-X-ray radia- 
tion at infrared wavelengths. Clumpy torus models pre- 
dict near ly isotropic emission in the MIR at wavelengths 
>12/xm (|Nenkova et all 120081 ). Thus, MIR-based surveys 
(or the combination of MIR and data at shorter wave- 



tion missed by hard X-ray surveys (e.g. 


Daddi at alj|2007l; 


Fiore at al. 2008; Gcorgantopoulos et al. 


20081: Fiore at alj 


20091: ISevergnini et alj 20121). For example, it has been 



claimed that objects showing excess emission at ~24^im 
over that expected from star formation, termed "infrared- 
excess galaxies" , might host heavily obscured and Co mpton- 
thick AGN (e.g. iFiore at al.ll200Sl ; iFiore at al.ll2009t ). How- 
ever the exact contribution of heavily obscured AGN to 
the infra red-excess galaxy popu lation remains an open is- 
sue (e.g. lAlexander et all 1201 j ). Several MIR-based AGN 
selection techniques have been developed with data from 
the Spitzer Space Telescope Infrared Array Camera (IRAC; 
I Fazio et~al I I2004J) using colours an d power-law selection 
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. 12004 IStern et al.j 120051 ; lAlonso-Herrero et al.1 
l200rj ; iDonlev et al.(l200Sl . I2012T ). These techniques are very 
effective and reliable. Galaxies dominated by AGN emis- 
sion typically exhibit a characteristic red power-law spec- 
tral ener gy distribution (SEP) in th e MIR (/„ oc v a with 
a ^-0.5; lAlonso-Herrero et al.ll200rj ) . Thus, MIR power-law 
selec tion provides the c leanest samples of luminous AGN 
(e.g. IDonlev et alj [20081 ). However, this technique is very 
sensi tive to the reliabilit y of the estimated photometric er- 
rors (|Donlev et al.ll2012l ). 

The Wide-field Infrared Survey Explorer (WISE) has 
now completed the first sensitive (~100-1000 x deeper than 
IRAS ) coverage of the entire sky in the MIBFI lWright et al.1 
2010). Several colour-based regions, aimed at identifying lu- 
minous AGN, have already been proposed. These works have 
shown that WISE can robustly separate AGN from normal 
galaxies and star s (e.g. lAssef et al.ll2010l ; |Jarrett et al.ll201ll ; 
IStern et al.ll2012l ). WISE will be extremely efficient in identi- 
fying the rare highly luminous AGN up to the crucial epoch 
when the accretion power of the Universe peaked (z~l-2). 
The all-sky WISE survey will complement the deep Spitzer 
surveys, aimed to characterize the accretion phenomenon in 
the distant Universe. 

This paper presents a highly reliable and complete MIR- 
based colour selection of AGN with WISE. Our technique 
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is designed to select objects with red MIR power-law SEDs 
and properly accounts for the estimated typical errors in 
the photometry and deviations of the MIR SEDs from a 
pure power-law. The AGN wedge is defined using the wide- 
angle Bright Ultra-hard XMM-Newton Survey (BUXS; Ma- 
teos et al. 2012c, in preparation). This survey is one of the 
largest complete flux-limited samples of bright "ultra-hard" 
(4.5-10 keV) X-ray selected AGN to date. Surveys such as 
BUXS are extremely efficient in selecting AGN bright enough 
for reliable optical identifications and for detailed studies of 
their propert ies and evolution (e.g. H BS28, ICaccianiga et al.1 
120041 ; HBSS. lDella Ceca et alfeoOSl ). BUXS covers the region 
of the AGN redshift-luminosity parameter space that WISE 
will sample. Thus, BUXS offers a unique opportunity to define 
a highly complete and reliable MIR-based AGN selection 
with WISE. Thanks to the optical spectroscopic identifica- 
tions available for ~97% of the BUXS objects, and the high 
quality X-ray spectra, we have maximized the completeness 
of our MIR selection without compromising its reliability. In 
a forthcoming paper we will present and discuss the main 
properties of the optical/near-IR/MIR SEDs of the AGN in 
BUXS (Mateos et al. 2012c, in preparation). 

This paper is organized as follows. Sections 2 and 3 
briefly summarize the data sets. In Section 4 we present our 
MIR selection of AGN candidates using the three shorter 
wavelength bands of WISE and the complete four bands, 
respectively and we discuss the completeness of the selection. 
We show the reliability of our AGN selection in Section 5. 
The results are summarized in Section 6. Throughout this 
paper errors are 90% confidence for a single parameter and 
we assume Q. m = 0.3, Q.a = 0.7 and Ho = 70kms~ 1 Mpc -1 . 



2 THE WISE INFRARED SURVEY 

WISE observed the entire sky in the MIR, achieving 5<r point 
source sensitivities better than 0.08, 0.11, 1, and 6 mjy at 
3.4, 4.6, 12, and 22 /im, respectively. The angular resolu- 
tion is 6.1", 6.4", 6.5", and 12.0" (FWHM), respectively, and 
the astrometric precision for high signal-to-noise (hereafter 
SNR) sources is better than 0.15" (Wright et al. 2010). We 
use here the March 2012 publicly available All-Sky Data Re- 
lease that covers >99% of the sky and incorporates the best 
available calibrations and data reduction algorithms (Cutri 
et al. 2012). 

In what follows we compute flux densities in the WISE 
bands using profile fitting photometry and the magnitude 
zero points of the Vega system: i ? „(iso)=309.124 Jy, 171.641 
Jy, 30.988 Jy, and 8.346 Jy for 3.4, 4.6, 12, and 22 ^tm, re- 
spectively. These values are computed with the flux correc- 
tion factors that correspond to a power-law sp ectrum (/„ oc 
u a ) w ith spectral index a=—l presented in IWright et al.l 
|2010l) . We note that using the flux correction factors that 
correspond to constant power-law spectra the difference in 
the computed flux densities would be less than 0.2% at 3.4, 
4.6, and 22 /xm and ~2% at 12 /im. We added a 1.5% uncer- 
tainty to the catalogued flux errors in all bands to account 
for the overall systematic uncertainty from the Vega spec- 
trum in the flux zeropoints. To account for the existing dis- 
crepancy between the red and blue calibrators used for the 
conversion from magnitudes to Janskys, an additional 10% 
uncertainty was added to the 12/im and 22/im fluxes (Wright 
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Figure 1. Redshift (top) and 2-10 keV luminosity (in log units, 
bottom) distributions of the BUXS type-1 and type-2 AGN. Lumi- 
nosities are rest-frame and are corrected for intrinsic absorption. 



et al. 2010). Throughout this paper we use monochromatic 
MIR flux densities (/„) in Janskys, unless otherwise speci- 
fied. 



3 THE BRIGHT ULTRA-HARD XMM-Newton 
SURVEY 

BUXS is one of the largest, amongst the existing XMM- 
Newton and Chandra surveys, complete flux-limited samples 
of bright (/4.5-iokev > 6 x 10 _14 erg s _1 cm -2 ) "ultra-hard" 
(4.5-10 keV) X-ray selected sources to date. BUXS is based 
on a subset of 381 high Galactic latitude > 20deg) 
observations from the se cond XMM-Newton serendipitous 
source catalogue (2XMM; W atson et all [20091 ). The sample 
is drawn from EPIC-pn observations with clean exposure 
times >10 ks and having good quality for serendipitous 
source detection (i.e. free of bright and/or extended X-ray 
sources). These observations were used to derive extragalac- 
tic source count distributions at intermediate fluxes, and 
there fore we have a good knowledge of the survey complete- 
ness (|Mateos et al.ll200Sl ). The total sky area of BUXS is 44.43 



deg 2 



The selection of sources in the 4.5-10 keV energy band 



was motivated by the need to reduce the strong bias against 
heavily absorbed AGN affecting surveys conducted at softer 
energies. The BUXS bright flux limit was intended to include 
only objects for which a reliable classification and redshift 
could be derived from optical spectroscopy. This also en- 
sures that we can derive accurate X-ray properties such as 
absorption and intrinsic luminosities from high-quality X- 
ray spectra (> few hundred counts). In this way we remove 
uncertainties associated with photometric redshifts and poor 
X-ray data quality (< a hundred counts). BUXS contains 258 
sources after removal of Galactic stars (<2%) and known 
BL Lacs (~1%). Optical spectroscopic identific ations have 
been obtained from the Sloan Digital Sky Survey (|Abazaiianl 
120091 ). the literature, and our ongoing follow-up campaign. 
At the time of writing, the spectroscopic identification com- 
pleteness is 97.3%. Of the 258 BUXS sources, 145 objects 
(56.2%) are identified as type-1 AGN (UV/optical emission 
line velocity widths ^1500 kms~ ) and 106 (41.1%) as type- 
2 AGN (UV/optical emission line velocity widths <1500 
kms~ or no emission lines). Seven sources (2.7%) remain 
unidentified. BUXS covers four decades in X-ray luminosity 
(~ 10 42 — 10 46 ergs~ 1 ), where the luminosities are computed 
in the 'standard' 2-10 keV rest-frame energy band and are 
corrected for intrinsic absorption. BUXS identifies sources out 
to z^2. Type-1 and type-2 AGN have mean 2-10 keV lumi- 
nosities of 1.5 x 10 44 ergs _1 and 2.2 x 10 43 erg s _1 , and mean 
redshifts of 0.7 and 0.3, respectively. Redshift and luminosity 
distributions are shown in Fig. [1] We note that BUXS sam- 
ples absorbed AGN in the Compton-thin regime (rest-frame 
intrinsic absorption Nh < 10 24 cm -2 ). 

To find the MIR counterparts of t he sources in BUXS w e 
used the cross-matching algorithm of IPineau et alj (|201ll ). 
The algorithm, which is based on the classical likelihood ra- 
tio, computes the probability of a spatial coincidence of the 
X-ray sources with their MIR candidate counterparts. MIR 
counterparts were found for 255 out of 258 (98.8%) sources 
(detection with SNR^5 in at least one of the WISE bands). 
The mean X-ray-MIR separation is <2 arcsec. The three 
objects without detection in the MIR are type-1 AGN with 
z~0.6-0.8. For one of these sources a blend of two WISE 
objects prevents us from identifying a unique MIR counter- 
part. For the other two sources, the expected MIR fluxes, 
from analysis of their optical/near-IR SEDs, suggest that 
they are too faint to be detected with WISE. 

In order to assess the completeness of our MIR selec- 
tion technique with respect to the overall AGN population 
we built a clean AGN parent sample by extracting all cata- 
logued WISE sources in the BUXS survey area, and we iden- 
tified the objects detected in the 2-10 keV ban d again using 
the cross-matching algorithm of IPineau" et al.1 (|201ll ). Here 
we assume that a detection in hard X-rays is a good tracer 
of unabsorbed and mildly ab sorbed AGN activity . We used 
the 2-10 keV source lists from lMateos et ail (|2008l ) that were 
derived with the source detection pipeline used in the 2XMM 
catalogue. The total number of objects in the BUXS area de- 
tected in the 2-10 keV X-ray band is 10265. The faintest 
sources have 2-10 keV fluxes of ~ 5 x 10 _15 erg s _1 cm~ 2 . 

We estimate a fraction of spurious matches of X-ray 
and MIR sources of <1% from the cross-matching of X-ray 
and WISE sources using a large offset in MIR coordinates 
(3 arcmin in either RA or dec). 
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Figure 2. Top: MIR colours for sources detected with SNR^5 at 3.4, 4.6, and 12 /im. Large symbols represent spcctroscopically identified 
BUXS AGN. Small cyan and grey symbols are WISE sources in the BUXS survey area with and without an X-ray detection at 2-10 keV, 
respectively. The 0<z<1.5 (Az=0.5) star-for ming tracks represen t M82 and the ULIRG Arp220. The 0<^^3 (Az=0.5) AGN tracks 
represent the infrared luminous AGN Mrk231 l|Polletta et al]|2008[) and a QSOl templ ate o btained by stitching together the 0.58-3.5Mm 
near-infrared and 3.5-24/im MIR quasar composite spectrum from lctlikman et al.l bOOrj ) and lHernan-Caballero &: Hatziminaogloul J201 lT) . 
Open symbols indicate z=0. Our AGN selection w edge and pow e r-law locus are the thick solid and dashed black lines, respectively. For 
comparison we show the AGN criteria defined by Ijarrett et al.l (12011!) and I Stern et al.l ||2012| ) , respectively (dotted-dashed purple and 
light blue lines). The error bars show the typical uncertainties in MIR colours at the SNR=5 limit (see Sec. 14. ll for details). Bottom: 
X-ray detection fraction of WISE objects across the colour-colour plane for bins containing at least 10 sources. Dotted-dashed (magenta) 
and dashed (yellow) contours indicate the density of WISE sources (normalized to the peak value) with and without X-ray detection, 
respectively. 

4 WISE SELECTION OF AGN CANDIDATES (3.4, 4.6, and 12/im in Sec. EQ 3.4, 4.6, 12, and 22^m in 

IN THE BUXS FIELDS Sec.g^J. 

To avoid spurious detections and objects with poorly con- 
strained photometry, in what follows we restrict ourselves 
to WISE sources detected with SNR^5 in all relevant bands 
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Table 1. Summary of the MIR selection of AGN candidates in 
the BUXS survey area. 



MIR wedge N WISE N WISE+X N wcdgo N wodgc+x 

(1) (2) (3) (4) (5) 

3- band 25206 1659 (6.6%) 2755 1062 (38.5%) 

4- band 2476 409 (16.5%) 516 245 (47.5%) 



galaxies, while the great majority of X-ray detected objects 
cluster near the power-law locus. 

Our MIR-based AGN selection technique is designed to 
identify objects with red MIR power-law SEDs. The AGN 
wedge is defined to include all objects with MIR colours 
expected for power-law SEDs with spectral index a ^—0.3, 
properly accounting for the typical errors in the photometry 
at the faintest MIR fluxes (error bars in Fig. [2j top). We 
then increase the size of the wedge towards red log(/4. 6 //3.4) 
colours (upper boundary) to include all BUXS AGN and the 
X-ray detected WISE objects throughout the BUXS survey 
area with such colours. In this way we account for deviations 
of the 3.4/jm to 12^m SEDs from a pure power-law. Our 
three-band AGN wedge is shown with the thick solid box in 
Fig. [3 The MIR power-law locus is defined by 

y = 0.315 x a; (1) 

where x = logio (|^) and y = logio (f^)- The top 
and bottom boundaries of the wedge are obtained by adding 
y-axis intercepts of +0.297 and —0.110, respectively. The 
MIR power-law a=— 0.3 bottom-left limit corresponds to 

y = -3.172 x x + 0.436 (2) 

These are the limits we have used throughout this pa- 
per. However, for most practical purposes using instead a 
bottom-left vertical limit, corresponding to s^0.120, pro- 
duces very similar results. 

In Vega magnitudes (in the WISE source catalogue 
magnitudes are reported in the Vega system), the AGN lo- 
cus is defined by 

y = 0.315 x x (3) 

where x'= [4.6] - [12] and y'= [3.4] - [4.6]. The top and bot- 
tom boundaries of the wedge are obtained by adding y-axis 
intercepts of +0.796 and —0.222, respectively. In this case, 
the MIR power-law a—— 0.3 bottom-left limit corresponds 
to 

y = -3.172 x x + 7.624 (4) 

The bottom-left vertical limit of the AGN wedge in magni- 
tudes corresponds to x-'^2.157. 

Our three-band AGN wedge identifies 2755 AGN can- 
didates in the BUXS area, of which 1062 (38.5%) are detected 
in X-rays (see Table[TJ). Out of the latter, 105 are associated 
with BUXS type-1 AGN and 38 with BUXS type-2 AGN (see 
Table [2]). We note that the X-ray detection fraction in the 
wedge increases with the depth of the X-ray observations as 
shown in Table(3] For example, in the BUXS area where the X- 
ray observations have exposures >40 ks, the X-ray detection 
fraction rises to 49.8%. For comparison, in the lMs CDF-S 
survey the X-ray detecti on fraction of IRAC power-la w AGN 
candidates was -50% {H onso- Herrero et all 120061 ) while 
this fracti on increased to ~8 5% in the deep er 2Ms CDF- 
N sur vey (|Donlev et alj|2007t ). Furthermore, iDonlev et al.l 
|2012l ) found that the X-ray detection fraction of IRAC MIR 
AGN candidates in COSMOS increased from 38% to 52% in 
the regions of deep Chandra coverage (X-ray exposures 50- 
160 ks). This is as expected, as long X-ray exposures are 
required to detect intrinsic ally less luminous and/or heav- 
ily obscured AGN (see e.g. iMateos et"al]|2005l ; iTozzi et afl 



Column 1: MIR AGN selection wedge identifier; Column 2: 
number of catalogued WISE sources in the BUXS survey area 
with significance of detection ^5 in the relevant bands; Column 
3: number (fraction) of WISE sources with an X-ray detection in 
the 2-10 keV band; Column 4: number of WISE sources in AGN 
wedge; Column 5: number (fraction) of WISE sources in AGN 
wedge and detected in X-rays. 



Table 2. Summary of the MIR selection of BUXS AGN. 



MIR wedge 


Opt. class 


NwiSE 


N wedge 


(1) 


(2) 


(3) 


(4) 




Type-1 


114 


105 


3-band 


Type-2 


81 


38 




No ID 


4 


3 




Typc-1 


63 


55 


4-band 


Type-2 


55 


21 




No ID 


2 


1 



Column 1: MIR AGN selection wedge identifier; Column 2: 
Optical class; Column 3: number of sources that are used to 
define the AGN wedge; Column 4: number of objects in AGN 
wedge. 



4.1 WISE three-band AGN wedge 

In the BUXS survey area there are 25206 sources detected 
with SNR^5 in the three shorter wavelength bands of WISE, 
of which 1659 have X-ray detections (see Table [TJ . Out of 
the latter, 114 are associated with BUXS type-1 AGN and 
81 with BUXS type-2 AGN (see Table 0). Fig. [5] shows the 
MIR log(/ 4 . 6 //3.4) vs. log(/i 2 // 4 .6) diagram for WISE ob- 
jects with and without an X-ray counterpart as small cyan 
and grey symbols, respectively. We also marked with large 
blue and red symbols spectroscopically classified type-1 and 
type-2 AGN in the BUXS survey. The dashed line illustrates 
the MIR power-law locus and the values for different spec- 
tral indices. Most BUXS objects, especially type-1 AGN, are 
clustered near the power-law locus, in a region in the MIR 
colour-colour plane well separated from the stellar locus 
(colours near zero magnitude) and the horizontal sequence of 
normal galaxies (lower right part of the diagram). Further- 
more, AGN have, on average, redder log(/4. 6 //3.4) colours 
than a pure power-law. This suggests s ome curvature in 
the observed 3.4/xm to 12/xm SEDs ( e.g. lElvis et al.1 1 19941 ; 
iRichards et~al]|2006l ; lAssef et al.ll201oi ). 

Fig. [2] (bottom) shows the fraction of WISE sources de- 
tected in X-rays across the colour-colour diagram and the 
distribution of objects with and without detection in X-rays 
(contours). There is a clear separation between these dis- 
tributions. The bulk of the MIR population not detected 
in X-rays overlaps with the horizontal sequence of normal 
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Figure 3. Equivalent to Fig. [2]for a MIR-bascd AGN selection using the complete four bands of WISE. Top: We include all objects with 
SNR^5 at 3.4, 4.6, and 12/xm and use different symbols depending on the 22/xm significance of detection. For SNR<2 at 22/xm we use 
la fluxes as upper limits. Symbols as in Fig. [2] 



l200fj ; IComastri et aljfeoill ; iBrightman et alj|2012h . Still, a 
substantial fraction of our MIR AGN candidates are unde- 
tected at 2-10 keV energies with the typical exposures in the 
2XMM catalogue. These objects that have the reddest overall 
log(/i2//4.6) colours in the AGN wedge, are the best can- 
didates to account for the most heavily obscure d /absorbed 
lumin ous AGN missed by hard X-ray surveys. lLacv et alj 
(|2007T ) presented the optical spectroscopic followup of a sam- 
ple of luminous AGN candidates selected on the basis of 
their IRAC MIR colours. They confirmed the AGN nature 
of 91% of the sources, with the majority of the objects be- 
ing identified as dust-reddened type-1 quasars and type-2 
AGN. Furthermore, a detailed study of the WISE MIR SEDs 



of [OIII]5007A-selected QS02s l|Reves et al.ll2008h strongly 
supports our hypothesis that many of the X-ray undetected 
sources in the wedge are heavily-obscured very-luminous 
AGN (Mateos et al. 2012b, in prep.). 

4.2 WISE four-band AGN wedge 

Due to the significantly shallower depth at 22^im compared 
with the first three bands, a selection that uses the 22/^m 
survey will be restricted by necessity to the brightest MIR 
objects. We have investigated, however, whether we can gain 
any additional information on AGN selection by using the 
complete four WISE bands. The number of MIR sources de- 
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Table 3. Dependence of the X-ray detection fraction of WISE 
sources in the 3-band AGN wedge as a function of the exposure 
time of the XMM-Newton observations. 



tcxp 


N wedge 


Nwedge+X 


(1) 


(2) 


(3) 


10-20 


1225 


400 (32.7%) 


20-30 


723 


298 (41.2%) 


30-40 


347 


135 (38.9%) 


40-50 


210 


110 (52.4%) 


>50 


250 


119 (47.6%) 


Total 


2755 


1062 (38.5%) 



Column 1: EPIC-pn exposure time interval of the X-ray 
observations in units of ks; Column 2: number of catalogued 
WISE sources in AGN wedge with significance of detection ^5 
in the three shorter wavelength bands of WISE. Column 3: 
number (fraction) of WISE sources in AGN wedge with an 
X-ray detection in the 2-10 keV band. 



tected with SNR^5 in all four WISE bands in the area of 
BUXS is 2476, of which 409 are detected in X-rays. Out of the 
latter, 63 are associated with BUXS type-1 AGN and 55 with 
BUXS type-2 AGN (see Tables Q] and [3J . Fig. [3] (top) shows 
the distribution of log(/4. 6 //3.4) vs. log(/ 2 2//i2) colours. The 
solid lines illustrate the four-band AGN selection wedge as 
in Sec. 4.1. In this case we increase the size of the wedge that 
would be required to account for the typical photometric er- 
rors (error bars in Fig. [3| towards both redder log(/ 2 2//i2) 
and log(/4.6//3.4) colours. In this way we increase the com- 
pleteness of the selection without compromising the reliabil- 
ity. The MIR power-law locus is defined by 



y = 0.50 x x 



where x 



(5) 

The top 



^/l2pm / y & \f3.iiJ, 

and bottom boundaries of the wedge are obtained by adding 
y-axis intercepts of +0.421 and —0.130, respectively. In this 
case we use a MIR power-law bottom-left limit of q=— 0.5 
that corresponds to 



-2.00 x x + 0.33 



(6) 



These are the limits we have used in Sec. 14.31 Using instead 
a bottom- left vertical limit, corresponding to a-^0.13, pro- 
duces very similar results. 

In Vega magnitudes the AGN locus is defined by 



V 



' = 0.50 x x 



{-) 



where x'= [12] - [22] and y'= [3.4] - [4.6]. The top and bot- 
tom boundaries of the wedge are obtained by adding y-axis 
intercepts of +0.979 and —0.405, respectively. In this case, 
the MIR power-law a=— 0.5 bottom-left limit corresponds 
to 



y = -2.00 x a;' + 4.33 



(8) 



The bottom-left vertical limit of the AGN wedge in magni- 
tudes corresponds to x '^1.76. 

Our four-band AGN wedge identifies 516 AGN candi- 
dates in the BUXS area detected with SNR^5 in all four 
WISE bands, of which 245 (47.5%) are detected in X-rays. 



BUXS type-1: 3-band wedge 
j»BUXS type-2: 3-band wedge 
^BUXS type-1: 4-band wedge 

BUXS type-2: 4-band we 
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Figure 4. Top: Fraction of BUXS AGN that meet our MIR selec- 
tion as a function of their intrinsic 2-10 keV luminosity (in log 
units). Filled symbols show selection using the three shorter wave- 
length bands of WISE and open symbols show selection using the 
complete four bands. The symbols indicate the mean luminosity 
of the sources in the bin. Bottom: Distributions of 2-10 keV lumi- 
nosity (in log units) for the BUXS typc-1 and type-2 AGN detected 
with SNR^5 in the three shorter wavelength bands of WISE. The 
luminosity (and redshift) distributions of the BUXS AGN detected 
in the three shorter wavelength bands of WISE and the complete 
four bands are consistent with each other. 



Out of the latter, 55 are associated with BUXS type-1 AGN 
and 21 with BUXS type-2 AGN (see Table HJ. 



4.3 AGN selection completeness 

Fig. HI (top) and TableHshow the fraction of BUXS AGN that 
meet our three-band MIR colour cuts as a function of their 
intrinsic 2-10 keV luminosity (filled symbols). The symbols 
indicate the mean luminosity of the sources in the bin. Rest- 
frame luminosities were derived from a detailed X-ray spec- 
troscopic analysis and are corrected for Galactic and any in- 
trinsic absorption. The completeness of our selection criteria 
is a strong function of luminosity for both type-1 and type- 
2 AGN. This result reflects the fact that objects with MIR 
colours not dominated by the thermal emission from the 
AGN will be missed by our selection. This effect is more im- 
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portant for low-luminosity AGN, especially if these sources 
are affected by large dust extinction at the shortest wave- 
lengths of WISE. In these objects the starlight from the host 
galaxy will dominate their MIR emission. Thus, it should be 
emphasized that the comparison of an MIR colour selection 
completeness for different classes of objects is only meaning- 
ful if the objects span the same range of luminosities. Taking 



Table 4. Dependence of the WISE 3-band AGN wedge complete- 
ness on the luminosity of the BUXS AGN. 



this into account, at L 2 



< 10 ergs 



■,% and 



39.1_ 9 5 % of the type-1 and type-2 AGN, respectively meet 
the selection. At £2-10 kev > 10 44 ergs _1 the MIR selection 
efficiency increases to 97. and 76.5+j|'4% for type-1 
and type-2 AGN, respectively. We estimated the most prob- 
able value for these fractions us ing a Bayesian approac h and 
the binomial distribution from IWall fc Jenkins! (|2008| ). The 
quoted errors are the narrowest interval that includes the 
mode and encompasses 90% of the probability (S. Andreon, 
private communication). It is important to note that the 
significantly smaller value of the selection completeness ob- 
tained for type-2 AGN at £2-10 kev < 10 44 ergs _1 compared 
to that for type-1 AGN is mainly due, as indicated above, 
to the different luminosity distributions of the two classes 
of AGN. The type-2 AGN population in BUXS is dominated 
by objects with £ 2 -iokev < 10 44 ergs _1 (~79% vs. ~39% for 
type-1 AGN; see Fig. [4] bottom). At such luminosities many 
AGN have relatively blue colours at the shortest WISE 
wavelengths (i.e. host-dominated) and thus lie outside the 
MIR AGN wedge. We expect this effect to be more impor- 
tant for type-2 AGN, as this class of objects is expected 
to show a higher degree of extinction at the shortest wave- 
lengths of WISE. For example, the clumpy torus models of 
iNenkova et al.l (|2008l ) predict nearly isotropic emission at 
wavelengths >12/^m. This could explain that even if we use 
the same luminosity range (as in Fig. [4| we obtain a se- 
lection completeness that is, within the uncertainties, still 
marginally lower for type-2 AGN than for type-1 AGN. 
This result still holds at luminosities L2-iokev > 10 44 ergs-\ 
where the relative contribution of the host galaxy to the 
MIR emission should be small. However, as the number of 
type-2 AGN in BUXS at such luminosities is small (16 ob- 
jects), the difference could be due in part to small number 
statistics. Thus, we conclude that our 3-band AGN wedge 
is highly complete for both X-ray selected luminous type-1 
and type-2 AGN. 

Fig. SI (top) and Table[5]show the fraction of BUXS AGN 
that meet our four-band colour cuts as a function of their 
intrinsic 2-10 keV luminosity (open symbols). The complete- 
ness of the four-band wedge is somewhat smaller but com- 
parable, within the uncertainties, to that achieved with the 
three-band selection for both type-1 and type-2 AGN. In- 
deed, 76 out of the 118 BUXS AGN detected in the four WISE 
bands meet the four-band selection, while this number in- 
creases to 88 if the three-band selection is used instead. All 
the 15 three-band selected BUXS AGN that miss the four- 
band selection have log(/ 2 2//i2) colours significantly redder 
than those expected for a pure power-law SED. These ob- 
jects lie outside the four-band AGN wedge, in the region 
of the colour-colour plane occupied by normal star-forming 
galaxies. They likely miss the four-band selection because 
their 22^im emission comes from both AGN activity and 
intense star formation. On the other hand, only three BUXS 
AGN from the four-band selection lie outside the three-band 
wedge. We note that we obtain the same result if we in- 
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Ntypo-2 

(4) 


/type — 2 

(5) 


[40-42] 
[42-43] 
[43-44] 
[44-45] 
[45-46] 


12 
33 
51 
18 


41 7+22.1 
-20.0 

96.1+1° 

100-11.4 


11 
17 
36 
16 
1 


0.0+ 174 
17 6+ 175 

1,0 -11.3 

61.1±™ 
75-0+ii 
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Column 1: X-ray luminosity range in units of ergs 
(logarithmic units, 2-10 keV in rest-frame and corrected for any 
intrinsic absorption); Column 2: number of BUXS type-1 AGN in 
luminosity bin; Column 3: fraction of BUXS type-1 AGN in the 
WISE 3-band AGN wedge; Column 4: number of BUXS type-2 
AGN in luminosity bin; Column 5: fraction of BUXS type-2 AGN 
in the WISE 3-band AGN wedge. 



Table 5. Dependence of the WISE 4-band AGN wedge complete- 
ness on the luminosity of the BUXS AGN. 
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63 
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Column 1: X-ray luminosity range in units of ergs 
(logarithmic units, 2-10 keV in rest-frame and corrected for any 
intrinsic absorption); Column 2: number of BUXS type-1 AGN in 
luminosity bin; Column 3: fraction of BUXS type-1 AGN in the 
WISE 4-band AGN wedge; Column 4: number of BUXS typc-2 
AGN in luminosity bin; Column 5: fraction of BUXS type-2 AGN 
in the WISE 4-band AGN wedge. 



elude all BUXS AGN with a lower significance of detection 
at 22/im (open symbols and arrows in Fig. 02 top). There- 
fore, by requiring 22/^m detections we are not biased against 
AGN with pure power-law SEDs (i.e. the faintest objects at 
22/^m). Thus, by including the 22/im WISE band to select 
AGN candidates neither the completeness nor the reliability 
of the selection improves. 

We have checked that using non-contiguous bands to 
define the WISE MIR colours does not increase the com- 
pleteness of a three-band or a four-band selection. 



5 RELIABILITY OF THE THREE-BAND AGN 
WEDGE 

In the following we restrict ourselves to the three-band AGN 
wedge which, as shown in the previous section, provides the 
most complete selection of AGN candidates in the BUXS fields 
with WISE. 
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Figure 5. Predicted z=0-2 WISE colours of AGN/galaxy composite SEDs for our three-b and AGN wedge. Th e AGN fraction is defined 
between 1 and 10 urn . The star- forming templates represent the UL IRG IRAS 22491 f square, iPolletta et al,l2008ri . th e starburst M82 (star , 
iPolletta et al,|[2008l) , a normal star-forming spiral galaxy (triangle, iDale at al.ll20o3 ) ■ and an elliptical galaxy (circle, [polletta et al, |2008). 
Large symbols mark each family of purely star-forming templates at z=0. The AGN template is the QSOl template of |Pollettaet*al.l 
(2008). Additional extinctions of Ay = — 2 and Ay = — 20 are applied to the star-forming and AGN components, respectively. Our 
MIR AGN selection wedge and power-law locus are the thick solid and dashed black lines, respectively. The power-law locus is defined 
from a=— 0.3. The colour tracks of purely star-forming galaxies would enter our selection wedge at z>1.3. At such redshifts however, 
most galaxies are too faint to be detected at the relatively shallow 12^tm flux density limits. 



5.1 Comparison with templates 

A known limitation of MIR selection techniques is the 
contamination from galaxies without AGN activity where 



Stern et al 



Stern et al 



.g. Lacv et alj 


2004 


iDonlev et al.1 


2012 



2012). To assess the reliability of our three- 



band AGN wedge we show in Fig. [5] the expected WISE 
colour s of AGN with varyi ng host-galaxy contributions fol- 
lowing [Donle^et^al] (|2012T ). The composite SEDs were con- 
structed using the library of IPolletta et al.l (|2008[) and the 
norma l star-forming spiral galaxy template from lpale at al.l 
(|2002h . The solid and dashed lines illustrate the AGN wedge 
along with the power-law locus, respectively. We also ap- 
plied additional extinctions of Av=0-2 and Av=0-20 to the 
star-forming and A GN components, respectively, using the 
iDraine et al.l (|2005f ) extinction curve. The colour tracks of 
pure star-forming galaxies would enter our selection wedge 
at z>1.3. At such redshifts, however, most galaxies are too 
faint to be detected at the relatively shallow WISE flux den- 
sity limits (|Wright et al.ll2010l ; Uarrett et alj|201lf ). Indeed, 
in the region of the AGN wedge where we would expect con- 
tamination from normal star-forming galaxies, there is no 
clear excess of MIR X-ray undetected sources and the X-ray 
detection fraction remains high (see Fig. [2j . Thus it seems 
that at the adopted SNR^5 limit our AGN selection suffers 



from minimal contamination from high redshift pure star- 
forming galaxies. It is interesting to note that an important 
fraction of the X-ray detected WISE sources has log(/i 2 //4.6) 
colours bluer than those expected for a pure AGN SED. The 
WISE colour tracks of composite galaxies suggest that this 
is most likely due to both AGN and their host galaxies con- 
tributing to the observed emission in the WISE bands. On 
the other hand, we find a sharp decrease in the X-ray detec- 
tion fraction of WISE objects at log(/i 2 // 4 .6) >0. 7-0.8. At 
such red MIR colours we expect many objects to be heav- 
ily obscured AGN. However, the colour tracks of pure AGN 
indicate that there is a strong dependence of the observed 
log(/i2//4.6) colour with redshift. This suggests that the pop- 
ulation of WISE objects with log(/ 12 // 4 .6) >0. 7-0.8 could be 
a mixture of heavily absorbed AGN and objects at high red- 
shifts (z>l-1.5). 



5.2 Comparison with other WISE selection 
techniques 

Within the BUXS survey area our three-band MIR colour 
selection identifies 2755 AGN candidates, of which 1062 
sources (38.5%) are detected in X-rays. For comparison, 
the X-ray detection fraction of WISE objects that meet the 
Ijarrett et al.l (|201ll ) colour-based AGN selection (indicated 
in Fig. [2]) is 33.4% (see Table 0. At red MIR colours their 
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Table 6. Comparison of our AGN wedge with other WISE selec- 
tion techniques. 



MIR wedge 


N wedge 


Nwedge+X 


(1) 


(2) 


(3) 


Mateos 3-band 


2755 


1062 (38.5%) 


Jarrett+11 


3301 


1102 (33.4%) 


Stern+12 


3946 


1254 (31.8%) 



Column 1: MIR AGN selection criteria; Column 2: number of 
catal ogued WISE AGN c andidates in the BUXS survey area. For 

the ljarrett et all (l201lh wedge we selected only MIR objects 
with significance of detection ^5 in the th ree shorter wavelen gth 
bands of WISE as in our analysis. The lStern et al. (2012) 
selection only requires an MIR detection at 4.6 /an, brighter 
than 160 fi3y (see Sec. 15.21 for details). Column 3: number 
(fraction) of WISE sources in AGN wedge with an X-ray 
detection in the 2-10 keV band. 

selection enters the sequence of low redshift normal galaxies 
incr easing the ex p ected number of contaminants. 

IStern et al.1 (|2012h proposed an AGN selection 
using a [3.4] -[4.6] colour cut ([3.4]-[4.6]>0.8 or 
log(/4. 6 //3.4) > 0.06) and a 4.6^im flux threshold of 
160^iJy. Their argument was that the inclusion of the longer 
wavelength WISE data would increase the reliability of the 
AGN selection but at the cost of reducing the completeness. 
Fig. [2] and Fig. [5] show tha t the [3.4]-[4.6] colour cut 
proposed by IStern et al.l (|2012l ) also enters the locus of low 
redshift normal galaxies at red log(/i 2 //4.6) colours. This 
could reduce the reliability of their selection as suggested 
by the lower X-ray detection fraction of WISE sources in 
the BUXS area that meet their criteria (31.8%; see Table [6]). 
On the other hand, using the dee p (60-160 k s ) Cha ndra 
data available in the COSMOS field. lStern et al.1 (|2012D find 
that 87% of of their WISE AGN candidates are detected 
at X-ray energies, suggesting minimal contamination from 
normal galaxies. 

The XMM-Newton pointings used to build BUXS span a 
broad range of ecliptic latitudes and thus the depth of the 
WISE survey varies across the BUXS fields. The lStern et al.l 
lj2012T ) AGN selection was defined using the COSMOS field 
located at low ecliptic latitude and thus, the MIR data are 
close to the minimum dept h of the WISE sur vey. At such 
shallow depths and using the lStern et~aH l|2012l ) 160/iJy flux 
threshold at 4.6^m, many of the star-forming contaminants 
will be too fai nt to be det e cted, thereby increasing the re- 
liability of the IStern et all (j2012T l selection. Although BUXS 
likely samples fainter objects than those best targeted at the 
shallow depth of WISE, we find that ~98% of the AGN have 
MIR detections with SNR^5 at 3.4/^m and 4.6/^m. This frac- 
tion only decreases to ~77% if we require 12/mi detections 
with SNR^5. In conclusion, over the range of MIR depths 
of the WISE survey in the BUXS fields, our proposed selec- 
tion suffers less contamination from star-forming galaxies 
than provided by a simple [3.4] — [4.6] color cut, while only 
marginally reducing completeness. 

We have investigated the impact of increasing the MIR 
significance of detection of the AGN candidates. If we re- 
quire 3.4/^m and 4.6^im detections with SNRJSlO, 44% of our 
MIR AGN candidates with log(/ 12 // 4 . 6 );?0.6 will be missed 
(21% of the objects with an X-ray detection), while this frac- 
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Figure 6. MIR colours of BUXS AGN as a function of intrinsic 
2-10 keV luminosity. Solid and dashed lines illustrate the AGN se- 
lection wedge and power-law locus, respectively. Most BUXS type-2 
AGN are objects with L 2 -iokev < 10 44 ergs _1 (~79% vs. ~39% 
for type-1 AGN; see Fig. |4}. At such luminosities many AGN 
have relatively blue colours at the shortest WISE wavelengths 
(i.e., host-dominated) and lie outside the AGN wedge. 

tion is only 0.7% at log(/i 2 //4.6)^0.6 (and zero for objects 
with an X-ray detection). Therefore, the effect of increasing 
the threshold in SNR of the detections is that the result- 
ing MIR selection becomes increasingly similar to a hard 
X-ray selection. By using SNR^5 WISE sources, we gain in 
the identification of AGN candidates with the reddest MIR 
colours, while we do not reduce the reliability. 

It is clear that our MIR colour selection is a good com- 
promise between completeness and the crucial high relia- 
bility required to obtain a clean sample of powerful AGN 
at the different depths of the WISE survey. Furthermore, 
going down to detections with SNR^5, we reach a much 
higher efficiency of detection of the AGN population in the 
reddest MIR colours, many of which could be heavily ob- 
scured/extincted AGN. 



5.3 Trends in the wedge with the AGN luminosity 

In the previous sections we presented the WISE expected 
colours for AGN/star- forming galaxies derived from a li- 
brary of templates. Here we investigate in more detail the 
colour trends in the three-band AGN wedge for the BUXS 
type-1 and type-2 AGN. Fig. [6] shows the MIR colours of 
the AGN in BUXS as a function of their X-ray luminosity. 
As expected, there is a strong dependence of log(/4. 6 //3.4) 
on the X-ray luminosity, with less luminous sources hav- 
ing the bluest MIR colours. For less powerful AGN the 
host gala xy can substantially contribute to the MIR emis - 
sion (e.g. iBuchanan et al 1 120061 : [Al onso-Herrero et al.ll2008T ). 
Host galaxy dilution is expected to be more important in 
type-2 AGN, especially at low luminosities, as type-2 AGN 
should show a higher degree of extinction at the shortest 
wavelengths of WISE. These objects, with MIR colours con- 
sistent with normal galaxies, lie outside of the AGN wedge 
and thus, they are missed from a pure MIR selection. Fig. [2] 
and Fig. [6] both show that most AGN in BUXS missed by 
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the 3-band wedge are type-2 AGN. As noted in Sec. 14.31 
in flux-limited X-ray surveys, such as BUXS, type-2 AGN 
are overall intrinsically less luminous than type-1 AGN (see 
Fig. [T] bottom). Thus, due to the strong dependence of 
the MIR selection completeness on the luminosity of the 
objects (see Fig. [31 top), our AGN wedge will preferen- 
tially pick out BUXS typc-1 AGN. However, at luminosities 
£2-10 kev > 10 44 ergs _1 , where the AGN is expected to dom- 
inate the MIR emission (unless it is heavily absorbed), the 
completeness of the selection of type-1 and type-2 AGN is 
comparable within the uncertainties. At such luminosities 
both type-1 and type-2 AGN are preferentially located above 
the power-law locus (see also Fig. [5]). This indicates that the 
observed 3.4/xm to 12/xm MIR SEDs of powerful AGN devi- 
ate from a pure power-law. 

We do not find a strong dependence of the log(/i 2 //4.6) 
colour on the X-ray luminosity. However, we find a large 
scatter in the distribution of log(/i2//4.6) colours, especially 
for objects at low luminosities, where the host galaxy sig- 
nificantly contributes to the MIR emission. A broad range 
of log(/i 2 //4.6) colours is expected for AGN with an impor- 
tant host galaxy contribution as the very wide 12/im filter of 
WISE is very sensitive to both prominent PAH emission fea- 
tures and silicate absorption (10/im) in star- forming galaxies 
over a broad range of redshifts. 



6 SUMMARY 

We present a MIR power-law based selection of luminous 
AGN candidates using the 3.4, 4.6, and 12 /im bands 
of the WISE survey. We defined an AGN wedge in the 
log(/4. 6 //3.4) vs. log(/i 2 //4.6) colour-colour diagram using 
the Bright Ultra-Hard XMM-Newton Survey (BUXS). This 
is one of the largest complete flux-limited samples of bright 
(/4.5-iokev > SxKr^crgs^cm- 2 ) "ultra-hard" (4.5-10 
keV) X-ray selected AGN to date. BUXS includes 258 objects 
detected over a total sky area of 44.43 deg 2 : 251 (97.3%) 
are spectroscopically identified and classified, with 145 
being type-1 AGN and 106 type-2 AGN. Our technique is 
based on a MIR power-law selection and properly accounts 
for the errors in the photometry and deviations of the MIR 
spectral energy distributions from a pure power-law. In 
flux-limited X-ray surveys, such as BUXS, type-2 AGN are 
intrinsically less luminous than type-1 AGN. Thus, due to 
the strong dependence of the MIR selection completeness 
on the luminosity of the objects, a MIR AGN wedge nec- 
essarily picks out BUXS type-1 AGN. However, at 2-10 keV 
luminosities above 10 44 ergs _1 the completeness of our MIR 
selection of type-1 and type-2 AGN is high and comparable 
for both types within the uncertainties. Our selection 
is highly complete at luminosities £2-10 kev > 10 44 ergs~ 
where our MIR wedge recovers ~97% and ~77% of the 
BUXS type-1 and type-2 AGN, respectively. We identify 
2755 AGN candidates in the 44.43 deg 2 BUXS survey area 
of which 38.5% have detection in X-rays. In the BUXS area 
where the X-ray observations have exposures >40ks, the 
X-ray detection fraction rises to 49.8%. This is reasonable, 
as long X-ray exposures are required to detect intrinsically 
less luminous and/or heavily obscured AGN. A substantial 
fraction of the MIR AGN candidates remain undetected 
at 2-10 keV energies with the typical exposures in the 



2XMM catalogue. These objects are the best candidates to 
account for the most heavily obscured/absorbed luminous 
AGN missed by hard X-ray surveys. Assuming that a 2-10 
keV X-ray detection is a good tracer of AGN activity we 
demonstrate that our WISE selection shows one of the 
highest reliability amongst those in the literature. This is 
crucial to obtain a clean MIR selection of powerful AGN. 
Furthermore, going down to a SNR^5 limit in the WISE 
flux densities, we substantially increase the efficiency of 
detection of AGN with the reddest MIR colours. We also 
investigate a WISE four-band AGN selection. We show, 
however, that by including the 22/xm WISE band neither the 
completeness nor the reliability of the selection improves. 
This is likely due to both the significantly shallower depth 
at 22/im compared with the first three bands of WISE and 
star-formation contributing to the 22/j.m emission at the 
WISE 22/xm sensitivity. 
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